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The ground state of the antiferromagnetic XY model witli a kagome lattice is known to be 
characterized by a well developed accidental degeneracy. As a consequence the phase transition in 
this system consists in unbinding of pairs of fractional vortices. Addition of the next-to-nearest 
neighbor (NNN) interaction leads to stabilization of the long-range order in chirality (or staggered 
chirality). We show that the phase transition, related with destruction of this long-range order, can 
happen as a separate phase transition below the temperature of the fractional vortex pairs unbinding 
only if the NNN coupling is extremely weak, and find how the temperature of this transition depends 
on coupling constants. We also demonstrate that the antiferromagnetic ordering of chiralities and, 
accordingly, the presence of the second phase transition are induced by the free energy of spin wave 
fluctuations even in absence of the NNN coupling. 

PACS numbers: 75.10.Hk, 64.60. Cn, 74.80.-g 



I. INTRODUCTION 

The antiferromagnetic XY model can be defined by 
the Hamiltonian 



H = Ji ^cos((^i - Lpj), 

NN 



(1) 



where Ji > is the coupling constant, ipi de- 
scribes the orientation of the classical planar spin 
Si = (cos(/?i, sin ipi), belonging to the site i of some regu- 
lar lattice, and the summuation is performed over the 
pairs of nearest neighbors (NN) on this lattice. The 
ground state of such model on a kagome lattice (Fig. 1) 
is known to have a huge accidental (that is not related 
to the symmetry of the Hamiltonian) degeneracy . 

For Ji > the minimum of the energy of each tri- 
angular plaquette is achieved when the three spins be- 
longing to it form the angles ±27r/3 with each other. In 
addition to the possibility of a simultaneous rotation of 
all three spins such arrangement is also characterized by 
the two-fold discrete degeneracy. When on going clock- 
wise around the plaquette the spins rotate clockwise (an- 
ticlockwise), the plaquette can be ascribed the positive 
(negative) chirality a = ±1. In the ground state of the 
antiferromagnetic XY model with triangular lattice the 
plaquettes with positive and negative chiralities regularly 
alternate with each other |^ . 

In any ground state on a kagome lattice the variables ipi 
analogously acquire only three values which differ from 
each other by 27r/3. However, the requirements for the 
arrangement of chiralities are less rigid than on triangu- 
lar lattice and accordingly the ground state in addition 
to the continuous U{1) degeneracy (related with an arbi- 
trary simultaneous rotation of all spins) is characterized 
by a well developed discrete degeneracy (of the same type 
as in the 3-state antiferromagnetic Potts model) leading 
to a finite residual entropy per site [|l|,^ . 

The accidental degeneracy persists if the interaction 
function in Eq. (Ill) differs from the pure cosine (remain- 



ing even), but is removed by the presence of interactions 
with more distant neighbors. In particular, for the ferro- 
magnetic sign ( J2 < 0) of the next-to-nearest neighbors 
(NNN) interaction the minimum of 



NN 



E 

NNN 



cos{ipi - ipj) (2) 



is achieved in one of the so-called Vs x y/3 states @] 
with a regular alternation of positive and negative chi- 
ralities. An example of such state is shown in Fig. 2(a). 
Here and further on we use the letters A, B and C to 
denote three values of fi which differ from each other by 
±27r/3 (for definiteness let us assume ^ (fiA + 27r/3, 
'PC = 'Pa+4:Tt/3). This state has the same structure as the 
ground state of a planar antiferromagnet with triangular 
lattice (or, to put it more accurately, can be obtained by 
the natural truncation of it). On the other hand the NNN 
interaction of the opposite sign ( J2 > 0) favors the fer- 
romagnetic arrangements of chiralities, which is achieved 
in the so-called q = states Q, see Fig. 2(b). 

For both signs of the NNN interaction the degeneracy 
of the ground state is reduced to U{l)x Z2 [1). This sug- 
gests the possibility of the two phase transitions, one of 
which is associated with unbinding of vortex pairs the 
Berezinskii-Kosterlitz-Thouless (BKT) transition [|-^] 
and the other with proliferation of the Ising type do- 
main walls. The number of the systems with the same 
degeneracy of the ground state includes, in particular, 
the antiferromagnetic XY model with triangular lattice 
[p|,p^ and the fully frustrated XY model with square 
lattice ||ll|,|l|. 

For the case of weak NNN interactin (IJ2I ^ Ji) the 
energy of a domain wall (per unit length) is proportional 
to I J2I and the logarithmical interaction of vortices to Ji. 
If domain walls and vortices were completely independent 
one would then expect the phase transition related with 
breaking of the discrete symmetry group to take place at 
lower temperature than the BKT transition. 
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In the present work we analyse the mutual influence 
between different classes of topological excitations in the 
antiferromagnetic XY model with a kagome lattice and 
weak NNN interaction and demostrate that the phase 
transition related with the domain walls proliferation (i) 
can happen as a separate phase transition only for ex- 
tremely weak NNN interaction, (ii) is not of the Ising 
type and (iii) the temperature of this transition is not 
proportional to the strength of the NNN interaction as 
one could naively expect. We also show that at very low 
temperatures the free energy of spin waves leads to sta- 
bilization of the antiferromagnetic ordering of chiralities 
even in absence of the NNN interaction. 

The results can be of interest in relation with possi- 
ble presence of weak easy-plane anisotropy in Heisen- 
berg kagome antiferromagnets [p3[|, which is indirectly 
confirmed by recent investigations of susceptibility ||lj] 
in (H30)Fe3(S04)2(OH)6 . The other class of physical 
systems, which allows for experimental realization of the 
considered model, consists of Josephson junction arrays 
1 15 1 and two-dimensional superconducting wire networks 
1 16 1 in perpendicular magnetic field. In such systems the 
role of ipi is played by the phase of the superconducting 
order parameter, and equivalence with the antiferromag- 
netic XY model is achieved when the magnitude of the 
magnetic field corresponds to half-integer number of su- 
perconducting flux quanta per each triangular plaquette. 
A superconducting array, which can be described by the 
same model in absence of the external magnetic field, can 
be constructed with the help of so-called 7r-junctions [p7|. 



II. ZERO TEMPERATURE: THE EQUIVALENT 
SOLID-ON-SOLID MODEL 



It has been already mentioned in the Introduction that 
the set of the ground states of the Hamiltonian (|l|) is 
equivalent (up to a simultaneous rotation of all spins) to 
that of the 3-state antiferromagnetic Potts model [|l],|j . 
Any triangular plaquette of a kagome lattice has to con- 
tain some permutation of the three values of (fi {(pA , <^b 
and fc)i which can be identified with the three states 
(a = 1, 2, 3) of the antiferromagnetic Potts model. 

The degeneracy of such set of ground states can be dis- 
cussed in terms of the zero-energy domain walls. If one 
consideres a \/3 x Vs ground state [Fig. 2(a)], it turnes 
out possible to construct the state with the same energy 
by permutation (for example) of the form </jb '^^^ V'c 
inside any closed loop formed by the sites with ipi = ipA, 
etc. Such closed loop [a simplest example of which is 
shown in Fig. 3(a)] can be considered as the zero-energy 
domain wall separating two different -\/3 x y/3 states with 
the opposite signs of staggered chirality. Any domain wall 
with zero energy is formed by elementary links which 
have to join each other at the angles of ±27r/3 [Fig. 3(b)]. 
Each link separates two triangular plaquettes with the 
same chirality, that is with the opposite signs of staggered 



chirality. The states with infinite (unclosed) domain walls 
are also possible. 

There exists ||l^ the exact mapping of the set of 
the ground states of the 3-state antiferromagnetic Potts 
model onto the states of the solid-on-solid (SOS) model 
in which the "height" variables u(r) are defined on the 
sites r situated at the centres of hexagonal plaquettes of a 
kagome lattice. These cites are shown in Fig. 1 by empty 
circles. They form the triangular lattice we shall denote 
T. Each site of the kagome lattice can be associated with 
some bond rr' of T and each variable ipi = (pA , <Pb , fc 
with the Potts variable a(rr') = a(r'r) defined on this 
bond. 

Since each triangular plaquette of a kagome lattice 
should always contain three different variables ipA, <^b 
and ipc, one can associate them with three basic vectors 
(ai + a2 + a3 = 0, see Fig. 4) of some auxilary trian- 
gular lattice, which we shall denote Ta {a — \aLa\ being 
its lattice constant). The height variables u(r), which ac- 
quire the values u G 7^ , can be then introduced following 
the rule 



u(r') = 



u(r) 
u(r) 



^Q(rr') 
^Q(rr') 



for r' 
for r' 



(3) 



where are the three basic vectors of T (ei +62 + 63 = 
0), as soon as the value of u(ro) is chosen for an arbitrary 
site ro [0. 

This defines the correspondence between the states of 
the antiferromagnetic Potts model and of the "vector" 
SOS model, in which the height variables u(r) € 7^ have 
to satisfy the constraint 



|u(r)-u(r')| = a 



(4) 



on all pairs of neighboring sites of T. By using the known 
properties of the exact solution [|l^ of the 3-state anti- 
ferromagnetic Potts model with external field coupled to 
staggered chirality Huse and Rutenberg [|l| have demon- 
strated (and also have confirmed this conclusion by nu- 
merical calculation) that such vector SOS model, in the 
partition function of which all allowed configurations of 
heights are counted with the same weight, is situated ex- 
actly at the point of the roughening transition, where (for 
In -r2| > 1) 

([u(ri) - n{v2)f) « ^ In |u(ri) - u(r2)|. (5) 

Therefore any additional perturbation supressing the 
fluctuations will lead to transition of the system into the 
flat phase, in which the fluctuations of u are convergent. 

According to constraint (^) the variables u(r) on neigh- 
boring sites have to be different from each other, so that 
even the most flat state is formed by the regular alterna- 
tion of three different values of u. The transition into the 
flat phase can be more transparently discussed in terms 
of the variables 



^(p^) ^ u(r) + u(rO + u(r'0 



(6) 
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describing the average height at each of the plaquettes of 
T. The variables n(R) are defined at the sites R of the 
honeycomb lattice which is dual to T, and acquire the 
values n(R) <E Tia, where Ha is the honeycomb lattice 
which is dual to 7^ (Fig. 4). In terms of the original spin 
variables the flat states [in which all variables n(R) arc 
equal to each other] correspond to \/3 x \/3 states, and 
zero-energy steps, the presence of which leads to their 
roughening, to the zero-energy domain walls separating 
different -s/S x V3 states. 

The large scale properties of the vector SOS model in- 
troduced above (and of its further generalizations) can 
be analysed with the help of the multi-component sine- 
Gordon model with the same symmetry. The (dimension- 
less) Hamiltonian of such sine-Gordon model can be cho- 
sen in the form 



H. 



SG 



[Vn(R)]^ 



y 



3 

E 

a=l 



cos[QQ,n(R) 



(7) 

The first term in Eq. (^ describes the effective stiffness 
(of entropic origin) which can be associated with the fluc- 
tuations of n(R), whereas the second term favors the val- 
ues of n(R) which belong to Tia- Here Qa are the three 
basic vectors of the triangular lattice reciprocal to Ta , so 

Qa = ^; Q1 + Q2 + Q3 



3a2 



0. 



Analogous Hamiltonian (with the opposite sign of the 
second term) and equivalent vector Goulomb gas have 
been investigated by Nelson in relation with dislocation 
mediated melting in two-dimensional crystals |2^. Al- 
ternatively the Hamiltonian of the form (|7|) can be in- 
terpreted as a simplified model for pinning of a two- 
dimensional crystal by a periodic substrate (cf. with 
Ref. [^). Note, however, that in contrast to real two- 
dimensional crystals, the accurate description of which 
requires to distinguish between compression and shear 
moduli, in our system the displacement n takes place 
in some auxilary space (and not in the real space) and, 
therefore, only one elastic modulus can be introduced. 

The renormalization group equations of Ref. , de- 
scribing the evolution of K and y with the change of the 
length scale L, in our notation can be rewritten as 

dK Stt o 



dl 

dy 

dl 



4ttK 



(8a) 
(8b) 



where I = In L. The corresponding flow diagram is 
schematically shown in Fig. 5, where Kc = l/(87r). It 
suggests that the roughening transition takes place when 
the renormalized value of the effective stiffness K is equal 
to Kc- The vector SOS model described above (which is 
known to be at the point of its roughening transition |^] ) 
can therefore be associated with some point belonging to 
the left separatrix. 



III. PHASE TRANSITION(S) ASSOCIATED 
WITH VORTEX PAIRS UNBINDING 



At finite temperature T other types of fluctuations (re- 
quiring finite energy) become possible, in particular for- 
mation of vortex pairs. Vortices are point-like topological 
excitations (the local minima of the Hamiltonian), the 
existence of which is related with the multivaluedness of 
the field ip. On going around each vortex ip experiences 
a continuous twist which adds to ±27r. At low temper- 
atures all vortices are bound in neutral pairs by their 
logarithmical interaction Q . 

With increase in temperature this interaction becomes 
renormalized due to mutual influence of different vor- 
tex pairs and becomes screened at the temperature Tbkt 
of the BKT transition (7|-H], which leads to dissociation 
of vortex pairs and exponential decay of correlations of 
exp («(/?) (in contrast to algebraic decay ^ at T < Tbkt)- 
The value of the helicity modulus F, describing the ef- 
fective stiffness of spin system with respect to infinitely 
small twist, at the temperature of vortex pairs dissocia- 
tion is known to satisfy the universal relation : 



Tbkt = ■^r(TBKT)- 



(9) 



Huse and Rutenberg |^] have argued that since at 
T — Q the antiferromagnetic XY model with kagome lat- 
tice is characterized by the long range order in exp(i3(/3) 
rather than in exp(i</5), the phase transition in this system 
should consist in unbinding of pairs of fractional vortices 
with topological charges ±1/3 and not of the ordinary 
(integer) vortices. The strength of the logarithmical in- 
teraction of fractional vortices is decreased by the fac- 
tor of 9 in comparison with that of integer vortices , 
therefore relation (0) should be replaced by 



(10) 



where Tpv is the temperature of the phase transition, 
associated with unbinding of pairs of fractional vortices. 

The value of F in any ground state of the antifer- 
omagnetic XY model with kagome lattice is equal to 
To = (V3/4)Ji. Substitution of Fo into Eq. ^ [instead 
of F(Tfv)] allows to obtain for Tpv the estimate (from 
above) of the form 



Ji « 0.075 Ji, 



(11) 



which turns out to be in reasonable agreement with the 
results of numerical simulations by Rzchowski |2^, who 
by using two different criteria has found Tpv ~ 0.070 and 
Tpv ~ 0.076. In Ref. [M the same estimate for Tpv has 
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been obtained in a less straightforward way with the help 
of the duality transformation |^,^ . 

The fractional vortices cannot exist by themselves (in 
absence of domain walls). A fractional vortex appears at 
every point where elementary links forming a domain wall 
meet each other at a wrong angle (tt/S or tt instead of 
27r/3). The same happens in the antiferromagnetic XY 
model with triangular lattice [Q, the ground state of 
which also has \/3 x %/3 structure. Fig. 6(a) shows an ex- 
ample of a domain wall containing one such special point. 
It separates the domain wall into two segments, one of 
which is formed by the sites with ipi — ipc and the other 
by the sites with ipi = ipB- 

When crossing the first segment the state to the right 
of the wall should be obtained from the state to the left 
by the permutation of A and B, whereas for the second 
segment the state to the right should be obtained by the 
permutation of A and C. This introduces the descrep- 
ancy of 27r/3 which can be localized on a semi- infinite 
line terminating at the special point (for example on the 
line X-Y-Z). In order to locally minimize the energy X, 
Y and Z should be replaced by A, B and C respectively 
when going from above, and by B, C and A when going 
from below. The misfit of 27r/3 has to be compensated 
by a continuous twist of (p, which is equivalent to the 
fractional vortex with the topological charge —1/3. 

In terms of the vector SOS model each fractional vor- 
tex corresponds to the point on going around which the 
height variable n changes by An with |An| = a. That 
means that at each fractional vortex a step with the 
height An (or, to put it more presicely, a set of steps 
with the total height An) terminates or begins. Accord- 
ingly the fluctuations of the SOS model provide the ad- 
ditional contribution to the interaction of the fractional 
vortices related with the difference in entropy between 
the configurations with different positions of step ter- 
mination points. At the point of roughening transition 
of the SOS model, as well as in the rough phase, this 
additional interaction (which can be expressed in terms 
of the correlation function of the dual XY model |27|) 
is also logarithmic. Its presence shifts Tpv upwards and 
diminishes the mutual influence between fractional and 
integer vortices. It is known [2^ that such mechanism in 
principle can lead to appearance of the separate phase 
transition, associated with unbinding of pairs of integer 
vortices, at temperatures above Tpv- 

On the other hand at finite temperatures the equiva- 
lence with the SOS model is no longer exact. One has to 
remember that the whole multitude of what we describe 
as flat states of the vector SOS model in terms of the 
original spin variables corresponds (for given (^a) to only 
six different \/3 x \/3 states, which can be obtained from 
the state shown in Fig. 2(a) by all possible permutations 
of A, B and C. In Fig. 4 the sites of TCa, which corre- 
spond to equivalent states in terms of (pi, are designated 
by the same numbers. In the zero-temperature partition 
function of the vector SOS model the properties of the 



zero-energy domain walls separting such states (in partic- 
ular, two closed loops formed by such walls cannot cross 
each other but can be situated inside each other or touch 
each other) allow to count them as different states of the 
SOS model At finite temperature it becomes pos- 

sible for a set of steps separating two physically equiva- 
lent states to terminate at the point where all these steps 
merge together |Q . The energy Ed of such defect is fi- 
nite and proportinal to Ji: Eo — cjjJi, where cjj is of 
the order of one. 

In terms of the multi-component sine-Gordon model 
(0) such defects correspond to dislocations of the field 
n, the (elementary) Burgers vectors of which (a = 
1, 2, 3) are given, as can be seen from Fig. 4, by 

bi=a3-a2, b2 = ai-a3, b3 = a2-ai. (12) 

An example of a dislocation is schematically shown in 
Fig. 6(b). It is formed by a neutral pair of fractional vor- 
tices that are sitting on two domain walls which cannot 
be transformed into a single domain wall. The letter X 
denotes the site on which ifi « ((^a + </'b)/2. In the vicin- 
ity of this site the values of (pi slightly deviate from those 
implied by the letters A, B and C. Successive application 
of the rule (||) along the perimeter of any closed loop 
surrounding point X sums up to An = a2 — ai . 

The renormalization of the dislocation fugacity z — 
expj—cpJi/T) with the change of the length scale can 
be described Q by 



dz 



X^z + 2nz^ 



where 



X,^2 = 2 - 4ttK 

An 



(13) 



(14) 



In the vicinity of the roughening transition [K « Kc) 
the exponent Az, describing the renormalization of z, is 
close to = 3/2, which corresponds to the fast growth 
of z. Comparison of Az with Ay = 2 — l/(47rif) shows 
that y and z are never simultaneously irrelevant. In that 
respect the situation is quite analogous to what is encoun- 
tered when considering the conventional (ferromagnetic) 
XY model with weak but relevant (low-order) anisotropy 

US- 

The presence of dislocations (or, to put it more accu- 
rately, of dislocation pairs) leads also to appearance in 
the right-hand side of Eq. (8a) of the additional (nega- 
tive) term proportional to z^. The presence of this term 
shifts the flow (see Fig. 5) from the separatrix to the area 
which corresponds to the rough phase of the SOS model. 

On the other hand the unrestricted growth of z under 
the renormalisation means that the system will contain 
the finite concentration of free dislocations, which trans- 
forms the rough phase of the SOS model into the disor- 
dered phase of the six-state model. The decay of correla- 
tions in this phase can be characterized by a finite corre- 
lation radius ^z, which can be found as the length-scale 
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at which (the renormaHzed value of z) becomes of the 
order of one. defines the scale at which the additional 
(entropic) interaction of the fractional vortices induced 
by the fluctuation of the domain walls is screened. The 
finiteness of closes even the hypothetical possibility for 
dissociation of pairs of integer vortices to take place as a 
separate phase transition at T > Tpv- 



IV. THE CASE OF THE FERROMAGNETIC 
NNN INTERACTION 

Inclusion into consideration of the interaction with 
more distant neighbors leads to removal of the acciden- 
tal degeneracy and stabilizes the states with either fer- 
romagnetic or antiferromagnetic ordering of chiralities of 
triangular plaquettes. 

In the case of the ferromagnetic NNN interaction 
( J2 < 0) the energy is minimal in one of the V3 X V3 
states with uniform staggered chirality [Fig. 2(a)]. Fig. 
3 shows two examples of a domain wall separating two 
different -^3 x states with opposite signs of stag- 
gered chirality. In the case of only NN interaction such 
domain wall (consisting of elementary links making an- 
gles of 27r/3 with each other) simply costs no energy. 
The presence of a weak ferrromagnetic NNN interaction 
makes the energy of such domain wall (per elementary 
link) i?DW equal to — 3J2 > 0. Here and further on we 
are interested only in the case | J2I Ji, when the values 
of the variables ipi remain close to those shown in Fig. 3 
not only away from the wall, but also in the vicinity of 
the wall. 

Note that such wall can fluctuate (make turnes, form 
kinks, etc.) without having to pay the energy propor- 
tional to Ji and, therefore, naively one could expect that 
the temperature of the phase transition, related with pro- 
liferation of such walls and leading to destruction of the 
long-range order in staggered chirality, should be deter- 
mined entirely by J2- Such conclusion, implicitely based 
on the comparison of the energy of an infinite domain wall 
with the negative entropic contribution to its free energy 
(the Peierls argument |Q), does not take into account 
that the presence of an infinite domain wall leads also to 
supression of the entropy (because it decreases the pos- 
sibilities for formation of closed domain wall loops) and 
in some cases does not work. 

In Sec. II we have discussed the properties of the vec- 
tor SOS model to which the antiferromagnetic XY-model 
with kagome lattice is equivalent at T = 0. In the parti- 
tion function of this model all allowed configurations of 
heights are counted with the same weight. In the case 
of the analogous model with a finite positive energy of 
a step (which corresponds to J2 < 0) the same partition 
function is reproduced in the limit of T — s- 00. For any 
small but finite ratio — J2/T > the SOS model is shifted 
from the point of roughening transition into the ordered 
(flat) phase. On the other hand we also know that at fi- 
nite temperatures the possibility of dislocation creartion 



tends to shift the same system into the disordered phase. 
One has to consider the competition of these two effects. 

In the vicinity of K = Kc the renormalization equa- 
tions (8) can be rewritten as 



dX 

m 

dl 



where 



X = 2\ l 



K 



Y 



1 

7ra ' 



1 

7!' 



(15a) 
(15b) 

(16) 



The solution of Eqs. (15) for arbitrary a allows one to 
find that the critical behavior of the correlation radius ^ 
in the vicinity of the transition is given |33] by 



In^ cx 



(ax) 



(17) 



where Aif is the deviation from the phase transition. In 
our problem for E'dw ^ T the ratio /S.K / is propor- 
tional to Ejy^ /T. 

The case of a = corresponds to Kosterlitz renormal- 
ization group equations [^ for the standard BKT tran- 
sition, which give V — 1/2. The case of a = +l/\/6 has 
been considered by Nelson [^^, who has found v = 2/5. 
The solution of the same equations for a = — 1/ gives 
17= 3/5. 

If the fugacity of dislocations z is so small that even 
when growing under renormalization it remains much 
smaller than one up to L ~ ^ [where the renormaliza- 
tion following Eqs. (15) stops anyway and fluctuations of 
n are frozen], the system remains in the ordered (flat) 
phase of the SOS model, that is in the phase with long 
range order in staggered chirality. On the other hand if 
Zi? (the renormalized value of z) manages to become of 
the order of one when the renormalized value of y is still 
small, the system finds itself in the disordered phase. The 
estimate for the temperature Tdw of the phase transi- 
tion separating these two regimes (and associated with 
the proliferation of domain walls) can be obtained from 
the relation z/?,(^) ^ 1, which is equivalent to 



ln(l/z)^A°lne 



In our problem 



ln(l/z) = cdJi/T 



and 



In^ 



T \ " 



DW 



(18) 



(19) 



(20) 



which leads to 
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Tdw ^ j Jl oc J2' Jl , (21) 

Away from the critical region the behaviour of 1^ can be 
described with the help of the self-consistent harmonic 
approximation , which gives V — Vq — 1. That means 
that with increase of J2/J1 the dependence ( pT| ) is re- 
placed by Tdw oc (J2Ji)"^/^. 

Note that the analysis which has led to Eq. ( pl| ) has 
been based on the assumption that all fractional vor- 
tices are bound in pairs, and, accordingly, is valid only 
for Tdw < Tpv- On the other hand the pairs of frac- 
tional vortices cannot dissociate at temperatures lower 
than Tdw, because for T < Tdw the fractional vortices 
in addition to their logarithmical interaction are bound 
also by domain walls (with a finite free energy per unit 
length) which connect them with each other. 

Therefore the two available possibilities are Tdw < 
Tpv and a single phase transition, whereas the scenario 
with Tdw > Jfv is impossible. Analogous conclusions 
have been earlier achieved in relation with hypothetical 
unbinding of fractional vortices in planar antiferromag- 
net with triangular lattice [Q. It is hardly surprising 
that the same conclusions are valid for the system, the 
ground state of which is practically identical to that of 
the antiferromagnetic XY model with triangular lattice, 
the only difference being that one quarter of sites is ab- 
sent. 

The proliferation of the low energy domain walls [of the 
type shown in Fig. 3(b)] leads to intermixing of six differ- 
ent states [which can be obtained from the state shown 
in Fig. 2(a) by arbitrary permutation of A, B, and C] and 
therefore should not be expected to be of the Ising type. 
Note that the domain walls are possible only between the 
states with the different signs of staggered chirality. The 
six-state model with analogous statistics of domain walls 
can be defined [ pO| by the partition function 

^6st= (n E ) n^(^i^-^i^') (22) 

V R tR=l/ NN 

where 

{1 for i = (mod 6) 

w fori =1,3, 5 (mode) (23) 
fori = 2, 4 (mode). 

The last line of Eq. ( p3| ) implies that the domain wall be- 
tween the states with the same parity of tji is impossible. 

Application of the duality transformation |35[| to the 
partition function (|2^) transforms it into analogous par- 
tition function with W(t) replaced by 

_ r 1 + 3w for i = (mod e) 

Vl^(t) = <^ 1 fori =1,2, 4, 5 (mode) (24) 

I 1 — 3w for i = 3 (mod e). 

The symmetry of W{t) corresponds to the so-called cu- 
bic model, the phase transition in the six-state version of 



which for W{1) > W{3) is known to be of the first order 
| p6| . The phase transition at Tdw in our system (at least 
when it happens at Tdw < Jfv) therefore also can be 
expected to be of the first order. 

Comparison of the estimate ( pl| ) with Eq. ( pT| ) shows 
that the fulfillment of the relation Tdw < Jfv requires 
< — J2 < Jmax, where Jmax can be estimated to be of 
the order of 10~^ Ji. For — J2 > Jmax there should be only 
one phase transition in the system, at which the prolifer- 
ation of domain walls is accompanied by the unbinding 
of all types of vortices. A detailed description of how it 
happens still remains to be constructed, but when the 
dissociation of pairs of fractional vortices is forced by the 
dissapearance of their linear interaction (mediated by the 
domain walls which connect them) at temperatures, for 
which their logarithmic interaction is already too weak, 
one can expect the value of the helicity modulus at Tc to 
be nonuniversal: 

2 ^(T^ 18 , , 

- < < - (25) 

Note that the estimate for Jmax has been found by taking 
the estimate (^l]) on its face value, that is without the 
unknown numerical coefficient, and therefore should be 
considered with great caution. 



V. THE CASE OF THE ANTIFERROMAGNETIC 
NNN INTERACTION 

For the antiferromagnetic sign ( J2 > 0) of the NNN in- 
teraction the minimum of the Hamiltonian (^) is achieved 
in one of the states with the ferromagnetic ordering of 
chiralites [Fig. 2(b)]. Such state also allows for construc- 
tion of a domain wall, the energy of which (per elemen- 
tary link) -Bd^v is proportional to the strength of second 
neighbor couphng: w 3J2, see Fig. 7(a). However, 

comparison of Fig. 7(a) with Fig. 7(b) shows that (in 
contrast to the case considered in Sec. IV) for J2 > 
the form of the state on the other side of the wall is 
uniquely defined by the orientation of the wall and is 
different for different orientations of the wall. The de- 
screpancy in tp that appears when crossing domain walls 
of different orientations should be taken care of by the 
fractional vortices which have to appear on all corners of 
domain walls (the same happens in the fully-frustrated 
XY model with square lattice [|l^). This makes impossi- 
ble the construction of a closed domain wall the energy of 
which is determined entirely by J2 and does not depend 
on Jl. 

For J2 ^ Jl a typical thermally excited defect (leading 
to the change of the sign of chirality) has the form of a 
long strip formed by two low energy domain walls [Fig. 
7(c)]. Like in Fig. e the letter X designates the sites with 
if « {ipA + Vb)/2. In the vicinity of these sites the values 
of other variables ipi slightly deviate from that shown in 



e 



the figure. Analogous strip defects are dominant at low 
temperatures in the frustrated XK-model with triangu- 
lar lattice and / = 1/4 or / = 1/3 

The energy of such defect is given by 2Eo 
where Eg = cqJi (cq ~ 0.55) is the energy of its termi- 
nation point and L its length. For J2 <C T <C Ji the 
average length (L) of such defects is given by the ratio 
T/{2E^^) 3> 1, whereas their concentration c is pro- 
portional to (L) exp(-2£'o/T). The relation c{L)'^ 1 
defines the temperature 



E. 



3Hn/E^^) 



(26) 



above which such defects no longer can be considered 
as independent. The same temperature can serve as the 
(lower bound) estimate for the temperature of the 
phase transition associated with proliferation of the do- 
main walls and leading to destruction of the long range 
order in chirality. For J2/J1 — > 



r» cx Ji/ln(Ji/J2) 



(27) 



Analogous estimate can be obtained by comparison 
of the domain wall energy E^^ with its entropy 
•^DW ~ 2exp(— £'x/r) due to possibihty of formation of 
kinks [Fig. 7(d)]. The energy of a kink Ek is very close 
to Eq. The requirement 
gives 



?+ 

DW 



Ek/H2T^w/Kw) 



(28) 



which is again the estimate from below. 

Like in the previous case (of the antiferromagnetic or- 
dering of chiralities) the proliferation of domain walls 
can take place as the independent phase transition only 
at temperatures lower than Tpv- Comparison of Eq. 
(p8|) with Eq. (pi]) shows that the fulfillment of relation 



^DW ^ ^FV requires J2 < ^maxj where Jj^j^x' '^^^ ^'^ ^^^i" 
mated as (10"^-;- 10~^) Ji. Also like in the previous case, 
the proliferation of domain walls is related with inter- 
mixing of six different states and therfore can hardly be 
expected to demonstrate the Ising type behavior. 



VI. SPIN WAVE FLUCTUATIONS 



in the second order in temperature Q . That is believed to 
be not sufficient for stabilization of a true long-range or- 
der related with chiralities. This conclusion does not take 
into account the percularities of the statistical mechanics 
of the considered system and has to be corrected. 

With the help of the numerical calculation (see Ap- 
pendix) we have found that the lowest order contribu- 
tion to the effective interaction of chiralities of neighbor- 
ing triangular plaquettes is of the antiferromagnetic sign 
(that is favors \/3 x \/3 state) and corresponds to 



E, 



(0) 

DW — ' 



Ji 



(30) 



where 7 « 2 • 10~^. Quantum and thermal anharmonic 
fluctuations in Heisenbcrg kagome antiferromagnet are 
also known to favour (at least locally) a planar state with 
%/3 X V3 structure @. The same can be said about the 
fluctuations of the order parameter amplitude in super- 
conducting wire networks | |4l[ |. 

Although E^^ defined by Eq. ( pO[ ) is always much 
smaller than the temperature and in the case of (for ex- 
ample) Ising model would be insufficient for appearance 
of the long-range order, in the considered system the sit- 
uation is qualitatively different. Substitution of Eq. ( |30| ) 
into Eq. (|l|) gives a finite value of Tdw induced by spin 
wave fiuctuations: 



(0) 



DW 



' 7' 



(31) 



which means that for 7^1 the long range order in 
staggered chirality would survive even up to T ~ Ji. 
However substitution of the numerically calculated value 
of 7 cited above produces an extremely low estimate: 



(0) 



DW 



10- Vi. 



Note that the ordering in staggered chirality is not- 
icable only at length-scales larger than the correlation 
radius ^. Substitution of Eq. ( pO| ) into Eq. ( [20| ) shows 
that for T < T^^ the behavior of ^(T) is given by 
In^ cx (Ji/7T)~. That means that at T -> there takes 
place a continuous reentrant phase transition into the 
phase without true long-range order in staggered chiral- 
ity. 



Another mechanism for removal of accidental degen- 
eracy (which is traditionally refered to as "ordering due 
to disorder" |Q) is related with the free energy of con- 
tinuous fluctuations (spin waves) Expansion of the 
Hamiltonian (|l|) up to the second order in deviations 
-tpi = ipi — (pi"-* of the variables ifi from their values ipf^"^ 
in some ground state gives the same answer 

NN 

for all possible ground states, which means that the dif- 
ference in the free energy between them can appear only 



VII. CONCLUSION 

This work has been devoted to investigation of the 
phase transitions in the antiferomagnctic XY model on a 
kagome lattice with the special emphasis on accurate con- 
sideration of mutual influence between different classes of 
topological excitations (fractional vortices and domain 
walls). In particular, we have shown that in the model 
with only NN interaction the additional interaction of 
fractional vortices related with the entropic contribution 
from zero-energy domain walls at flnite temperatures be- 
comes short-ranged. Therefore it can not interfere with 
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the BKT dissociation of fractional vortex pairs proposed 
in Refs. §] and 

For the case of a finite NNN coupUng (leading to 
removal of the accidental degeneracy) we have demon- 
strated that the phase transition related with prolifera- 
tion of the domain walls can happen as a separate phase 
transition below Tpv only for very weak NNN interaction, 
and have found how the temperature of this transition 
depends on Ji and J2 ■ These dependences are essentially 
different for different signs of the NNN coupling. The 
same results are also applicable for other mechanisms of 
removal of the accidental degeneracy, which lead to a fi- 
nite i?DW- Note that our analysis has been restricted to 
the case | J2 1 ^ Ji , so we have not considered the possibil- 
ity of the domain wall proliferation happening above the 
temperature of the ordinary BKT transition, associated 
with appearance of free integer vortices (like it happens 
in the case of triangular lattice | |lO| , ^ , ^ ) . 

Our conclusions are compatible with the results of 
the numerical simulations of Geht and Bondarenko 
who have found (for not too weak NNN interaction, 
I J2I > O.lJi) that the disordering of all degrees of free- 
dom in the antiferromagnetic XY model with kagome 
lattice takes place at the same temperature, the sin- 
gularities of the thermodynamic quantities being of the 
Ising type. Recently it has been shown |^ (for the case 
of triangular lattice) that when the domain wall pro- 
liferation happens as a continuous phase transition [at 
T/T{T) > 2rFv/r(rFv)], the dissociation of pairs of in- 
teger vortices has to take place at T < Tdw Since the 
same arguments are also applicable for a kagome antifer- 
romagnet with J2 < 0, it may be of interest to check the 
results of Ref. |^ with better accuracy. 

The long-range order in staggered chirality is favored 
also by the spin-wave fluctuations. Our analysis suggests 
that the antiferromagnetic XY model with kagome lat- 
tice and only NN interaction presents a unique exam- 
ple of a model without free parameter in which one of 
the phase transitions can be expected to happen at di- 
mensionless temperature of the order of 10~^. Therefore 
one can conclude that the numerical simulations of Rz- 
chowskii have demonstrated no evidence for selection 
of a single ground state down to T/Ji ^ 10~^ not be- 
cause i^DW oc T^/Ji is not sufficient for that, but simply 
because the temperature was not low enough. Compari- 
son with Eq. (^^ shows that if the effective interaction 
of chiralities, induced by the free energy of spin waves, 
would be of the opposite sign, the long range order in 
chirality would persist up to much higher temperatures. 

Experimentally the phase transitions discussed in this 
work can be observed in superconducting wire networks 
or Josephson junction arrays in the external magnetic 
field providing one-half of the superconducting flux quan- 
tum (j)o = hc/2e per each triangular plaquette of a 
kagome lattice. In such systems the removal of the ac- 
cidental degeneracy is related with the magnetic interac- 
tion of the currents and a finite width of the wires . 
The former of these mechanisms favors the ferromagnetic 



ordering of chiralities, whereas for the latter the effect de- 
pends on the width of the wires. 

Recent experimental investigation of the aluminum 
wire network with kagome structure [ p^ has demon- 
strated for (/) = 0o/2 the presence on the current- volatage 
curve of the regions corresponding to different mecha- 
nisms of dissipation, one of which (with an algebraic be- 
havior) can be associated with unbinding of vortex pairs 
and the other with spreading of domains with inversed 
chiralities |Q . The authors of Ref. have interpreted 
this as an evidence for the presence of two phase transi- 
tions. 
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APPENDIX: 

The lowest order contribution to the interaction of chi- 
ralities of neighboring triangular plaquettes induced by 
the spin wave free energy appears when the partition 
function of the Hamiltonian (|l]) is expanded up to the 
second order in 
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(Al) 



NN 



and then is averaged with the help of H^^\ The fourth- 
order term is the same for all the ground states and there- 
fore of no importance. 

The parameter describing the effective interaction 
of chiralities of neighboring triangular plaquettes (a and 
b): 

E{(Ja,(Jb) = K„(Ja(Tb (A2) 

can be then found by calculating the average of 



V = - 



T 



(A3) 



where 



27T 



= ^^Mii'i - ^0)' + (^-2 - V^i)' + (^0 - V'2)'] 

6 

(A4) 
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and exp ression for i/f ^ can be obtained by replacing in 
Eq. (A4) ipi by -03 and '02 by tp4. The indices from to 4 
are used here to denote the five sites belonging to a pair 
of neighboring triangular plaquettes as shown in Fig. 8. 

With the help of the Wick's theorem and symmetry 
arguments the average of V can be reduced to the form 

{V)^^{Gl-2GlGs + 2G,Gl-Gl) (A5) 



where 



and 



^3 — 9oi — 2^13; — 9oi ~ 2^14 



9ij = {{'P^ - "Pjf) 



(A6) 



(A7) 



describes the amplitude of flutuations of Lpi — ipj calcu- 
lated with the help of the harmonic Hamiltonian (p9|). 

The value of gij for the nearest neighbors (501) can be 
calculated exactly: 



301 



T 



(A8) 



whereas numerical calculation of the integrals over Bril- 
louin zone defining 1713 and 314 gives 



313 



3 ^ 



T 



(A9) 



where d w 0.0213. 

Substitution of Eqs. (^-(^) into Eqs. (|A5|)-([A^) 
then gives 



K„ = 



111 
2Ji 



where 



32 



5{1 + 12(5^) w 2.01 • 10" 



(AlO) 



(All) 



which leads to Eq. {Bu 
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FIG. 1. A kagome lattice (shown by filled circles) can be 
constructed by the regular elimination of one quarter of sites 
from a triangular lattice. It consists of triangular and hexag- 
onal plaquettes. 



FIG. 2. The structure of the ground states selected in pres- 
ence of interactions with further neighbors. The letters A, B 
and C correspond to the three values of ifii, which differ from 
each other by ±27r/3. (a) a -\/3 x y/3 state; (b) a q = state. 



FIG. 3. Two examples of a domain wall separating different 
states. 



FIG. 4. The triangular lattice Ta and its three basic vectors 
a^. The sites of the dual lattice Tla are shown by the num- 
bers from 1 to 6. The same numbers correspond to physically 
equivalent states. 



FIG. 5. The schematic flow diagram for Eqs. (8). The sys- 
tem with only NN interaction and T = can be associated 
with some point (shown by black dot) on the left separatrix. 
Dashed arrow shows how the flow is changed at T > when 
the contribution from z has to be taken into account. 



FIG. 6. (a) An example of a fractional vortex, (b) An ex- 
ample of a dislocation formed by a neutral pair of fractional 
vortices. 



FIG. 7. (a) and (b) two examples of a domain wall sepa- 
rating two q = states; (c) a typical finite size defect on the 
background of q = state; (d) a kink on a domain wall. 



FIG. 8. The numbering of sites used in the expressions for 



fff' and H, 



(3) 
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